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Abstract 

We study the structure of QCD at very large baryon density for an arbitrary number of flavors 
Nf. We provide evidence that for any number of flavors larger than Nf = 2 chiral symmetry 
remains broken at asymptotically large chemical potential. For N c = Nf = 3, chiral symmetry 
breaking follows the standard pattern SU(3)l x SU — > SU(3), but for Nf > 3 unusual patterns 
emerge. We study the case Nf = 3 in more detail and calculate the magnitude of the chiral order 
parameters (t/>^>) and (('tpip) 2 ) in perturbative QCD. We show that, asymptotically, (f/^) 1//3 is 
much smaller than ((^'0) 2 ) 1 ^ 6 • The result can be understood in terms of an approximate discrete 
symmetry. 



I. INTRODUCTION 



The phase structure of matter at non-zero baryon density has recently attracted a great 
deal of interest. In particular, it has been emphasized that quark matter at very high density 
is expected to behave color superconductor . The behavior of hadronic matter in 
this regime is of interest in understanding the structure of compact astrophysical objects 
and the physics of heavy ion collisions in the regime of maximum baryon density. Moreover, 
it was realized that matter at very high density exhibits many non-perturbative phenomena, 
such as a mass gap and chiral symmetry breaking, in a regime where the coupling is weak 
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. This means that some of the "hard" problems of QCD can be studied in a systematic 
fashion. 

At very high density the natural degrees of freedom are quasiparticles and holes in the 
vicinity of the Fermi surface. Since the Fermi momentum is large, asymptotic freedom 
implies that the interaction between quasiparticles is weak. In QCD, because of the presence 
of unscreened long range gauge forces, this is not quite true. Nevertheless, we believe that 
this fact does not essentially modify the argument However, as we know from the theory 
of superconductivity the Fermi surface is unstable in the presence of even an arbitrarily 
weak attractive interaction. At very large density, the attraction is provided by one-gluon 
exchange between quarks in a color anti-symmetric 3 state. QCD at high density is therefore 



expected to behave cLS cL color superconductor 

Color superconductivity is characterized by the breakdown of color gauge invariance. As 
usual, this statement has to be interpreted with care. Local gauge invariance cannot really 
be broken j^. Nevertheless, spontaneously broken gauge invariance is a useful concept. We 
can fix the gauge, introduce a gauge non-invariant order parameter and study its effect on 
gauge invariant correlation functions. The most important gauge invariant consequence of 
superconductivity is the appearance of a mass gap, through the Meissner-Anderson-Higgs 
phenomenon. The formation of a mass gap is of course also characteristic of a confined phase, 
ndeed, it is known that in general Higgs and confined phases are continuously connected 

Color superconductivity may also lead to the spontaneous breaking of global symmetries. 
It is this phenomenon that we wish to study in more detail. Broken global symmetries lead 
to the appearance of Goldstone bosons, and determine the low energy effective description 
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of the system. It is one of the remarkable properties of the color superconducting phase in 
Nf = 3 QCD that the pattern of broken global symmetries exactly matches that of QCD at 
low density. In this work, we would like to establish the global symmetry breaking pattern in 
QCD with Nf = 2, 3, . . . , 6. In particular, we wish to show that the result in the case Nf = 3 
is generic in the sense that for any number of flavors larger than 2, chiral symmetry is broken 
but a vector-like flavor symmetry remains. We also wish to study the most interesting case, 
QCD with Nf = 3, in more detail. We calculate the magnitude of the totally symmetric 
and anti-symmetric order parameters, and the magnitude of the chiral condensates {^ip) 
and {(ipip) 2 )- We comment on a number of unusual features of chiral symmetry breaking in 
QCD at large density. 



II. THE EFFECTIVE POTENTIAL 



In this section we wish to introduce a simple energy functional that captures the essential 
dynamics of QCD at high baryon density. We will use this functional in the following section 



in order to analyze the ground state of QCD with an arbitrary number o 



' flavors. 



12 1 . In this work 



Our construction is guided by the renormalization group analysis of 
we classified possible instabilities of the Fermi surface, and assessed their relative impor- 
tance, for small and short-range, but otherwise arbitrary couplings near the Fermi surface. 
It was found that the dominant instability corresponds to scalar diquark condensation. The 
analysis does not fix the color and flavor channel of this instability uniquely, because there 
are two equally enhanced interactions. One gluon exchange, which dominates for weak 
coupling, is attractive in the color anti- symmetric 3 channel, and favors one of these inter- 
actions. The same is also true for instanton induced interactions, that are likely to play an 
important role at moderate densities. As the couplings evolve towards the Fermi surface, 
the attractive interaction in the color 3 channel will grow, while the repulsive interaction in 
the color symmetric 6 channel is suppressed. The dominant coupling is a color and flavor 
anti-symmetric interaction of the form 

C = G(5 ac 5 bd - 5 ad 5 hc )(5 lk 5 3l - 5 u 5 jk ) {(^C^}) (fiC-yetf) - (C 75 <- C)} , (1) 

where a,b, . . . are color indices and i,j--- are flavor indices. In the following, we shall use 
the notation fC^f for the color-flavor structure of the interaction. 
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In full QCD, there are unscreened magnetic gluon exchanges and the interaction between 
quarks is not short range. This problem was first studied in [5]. We expect that this effect 
does not modify our results for the structure of the ground state, but only the magnitude 
of the gap. We check this explicitly in the case of Nf = 3 in section II VI In this section, 
we also study the effect of color-symmetric interactions. It was also pointed out that for 
asymptotically large chemical potential, the gap equation in higher partial wave channels 
becomes degenerate with the s-wave gap equation 0|. Subleading effects are expected to 
lift this degeneracy. Using the methods developed in ^| we have checked that the s-wave 
gap is bigger than higher partial wave gaps. We therefore expect the ground state to be an 
s-wave superconductor. The only exception is QCD with only one flavor, since in this case a 
color 3 order parameter cannot have spin zero. This situation is of physical interest for the 
behavior of real QCD with two light and one intermediate mass flavor. For m s lar ger than 



some critical value, QCD has a phase with separate pairing in the ud and s sectors jl3l \14\. 
In this work, we will not consider the phase structure of Nf = 1 QCD. 

In order to determine the structure of the ground state we have to calculate the grand 
canonical potential of the system for different trial states. Since the interaction is attractive 
in s-wave states, it seems clear that the dominant order parameter is an s-wave, too. We 
then only have to study the color-flavor structure of the primary condensate. We assume 
that the condensate takes the form 

W<?75^> = € ? • (2) 
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is a Nf x Nf matrix in flavor space and a N c x N c matrix in color space. The Pauli 



principle requires that is symmetric under the combined exchange (ai) <->• (bj). 

We calculate the effective potential using the bosonization method. For this purpose, 
we introduce collective fields Af- and A^ with the symmetries of the order parameter (j2J). 
We add to the fermionic action a term (4G) _1 /C^f A^A-y and integrate over the dummy 
variables A" b and A" b . We then shift the integration variables to eliminate the interaction 
term (J!}. So far, no approximations have been made. We now assume that the collective 
fields are slowly varying, and that A^- can be replaced by its vacuum expectation value. In 
this case, we can perform the integration over the fermionic fields and determine the grand 
canonical potential as a function of the gap matrix A"j\ 

The integration over the fermions is performed using the Nambu-Gorkov formalism. We 
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introduce a two component field ^ = ). The inverse quark propagator takes the form 



1 ' ii + iji-m K-A 



K ■ A (fi-i/L + m) T 
The grand canonical potential is now given by 



(3) 



V 



1. 
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fi(A) = -Tr log^S) + —A -K-A. (4) 
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In order to evaluate the logarithm, we have to diagonalize the mass matrix Ai = /C • A. Let 
us denote the corresponding eigenvalues by A p (p — 1, . . . , N c Nf). These are the physical 
gaps of the NfN c fermion species. If we neglect the quark masses m, the grand potential is 



(5) 



The momentum integral has an ultraviolet divergence. This integral can be regularized by 
introducing a cutoff A or, more generally, by including a form factor F(p). In this case it 
would seem that the properties of the grand potential depend on a number of parameters, 
such as the coupling constant G, the chemical potential /i and the cutoff A. In the weak 
coupling limit, however, the grand potential should only depend on the value of the gap on 
the Fermi surface, and not on the exact momentum dependence of the interaction. This 
can be made manifest by introducing a renormalized grand potential. Following the work 
of Weinberg jl^ , we have 



(6) 



Here, £ is a renormalization scale. The grand potential is independent of £, since the scale 
dependence of the first term is canceled by the scale dependence of the coupling constant 
G. The coupling constant satisfies the renormalization group equation discussed in jlll . l^ . 

The grand potential (jBJ) depends on N C (N C — l)Nf(Nf — l)/4 parameters. We minimize 
this function numerically. After we determine the matrix A° b that minimizes the grand po- 
tential we study the corresponding symmetry breaking pattern. Without superconductivity, 
there are Nj — 1 global flavor symmetries for both left and right handed fermions, as well 
as — 1 local gauge symmetries. Superconductivity reduces the amount of symmetry. 
In order to find the residual symmetry group we study the second variation of the order 
parameter 5 2 A/(58i58j), where 9i (i = 1, . . . ,Nj + A^? — 2) parameterizes the flavor and 
color transformations. Zero eigenvalues of this matrix correspond to unbroken color-flavor 
symmetries. 



III. COLOR SUPERCONDUCTIVITY IN QCD WITH iV c = 3 COLORS AND N j 
FLAVORS 



The results of our numerical study for QCD with three colors and Nf = 2, . . . , 6 flavors 
are summarized in Table 1. In the following, we will study each of these cases in more detail. 
The two flavor case is special. In this case, the order parameter 

only breaks color SU(3) — > SU(2). The chiral SU(2)l x SU(2)r symmetry remains unbro- 
ken. At this level, the Fermi surfaces of the up and down quarks of the third color remain 
intact. A careful study of the quantum numbers of the low ener gy s tates shows that chiral 
symmetry is realized in terms of massless protons and neutrons [13]. The proton and neu- 
tron are composites of the quark and Higgs fields. Subleading interactions can generate a 
gap for these states. The exact nature of this gap is hard to determine, even in the limit of 
very large chemical potential. 

For Nf larger than 2, the gauge symmetry is always completely broken, and all quarks 
acquire a gap. In addition to that, we find that the preferred order parameter always involves 
a coupling between the color and flavor degrees of freedom. This means that the original 
flavor symmetry is broken, but some vector-like symmetry which is a combination of the 
original flavor and color symmetries remains. 

In the case of three flavors we find that the preferred order parameter is of the form 

A% = A(5:6 b j -5?6 b i ). (8) 

This is the color-flavor locked phase suggested in Both color and flavor symmetry 
are completely broken. There are eight combinations of color and flavor symmetries that 
generate unbroken global symmetries. The symmetry breaking pattern is 

SU(3) L x SU(3) R x U(l) v -> SU(3) V . (9) 

This is exactly the same symmetry breaking that QCD exhibits at low density. This has 
led to the idea that in QCD with three flavors, the low and high density phases might be 
continuously connected j^]. We also note that the quark mass gaps fall into representations 
([8] + [l]) of the unbroken vector symmetry. Note that in the present analysis, which only 
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takes into account the leading interaction, the order parameter is completely anti-symmetric 
in both color and flavor. In the case Nf = N c , however, there is a more general order 
parameter Afj* = A 1 5"5j+A 2 5"(5 l fe which has both symmetric and anti-symmetric components 
and has the same residual symmetry. We will discuss this situation in more detail in the 
next section. 

In the case Nf = 4 the numerical results indicate that the most favorable order parameter 
is given by 

Ag> = Ae o6c ^. = Ae abc (e l]C + 6 ic 6, 4 - 5 jc 5 H ), (10) 

where rj^ is the 't Hooft symbol. There is a second, degenerate, solution where we replace 
Vij ~^ Vij = e ija — fiiafiji + Sj a 8i4. The 't Hooft symbol describes the isomorphism SU(2) x 
ST/ (2) = 0(4). Consider the 0(4) generators (M p<7 )^ = i(5ffi-5ffi). The 0(4) generators 
can be decomposed into SU(2) L xSU(2) R generators ^^(M^) pa and 77° (M**")^. Using this 
result one can show that the Nf = 4 order parameter (fTUj) realizes the symmetry breaking 
pattern 

SU(4) L x SU{A) R x U(l) v -> SU(2) V x SU(2) V . (11) 
The generators of the SU(2) x SU(2) symmetry are 

(e ahc N bc + \r)%M^ , (W fec + \% V M^ . (12) 

Here, N bc are 0(3) C SU(3) color generators and M pu are 0(4) C St/ (4) flavor generators. 
We note that for Nf = 4, and indeed for any Nf larger than three, we cannot realize the 
same symmetry breaking that we have at low density. This means that even though chiral 
symmetry remains broken at very large density, there has to be a phase transition that 
separates the high and low density phases. 

In the case Nf = 5 we were unable to find a compact expression for the energetically 
favored order parameter. The numerical results indicate that Nf = 5 QCD realizes the 
symmetry breaking pattern 

St/(5) L x SU(5) R x U(l) v -> SU{2) V . (13) 

The residual symmetry is smaller than for any other number of flavors. This is also reflected 
in the fact that the condensation energy is comparably small. 
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If Nf is a multiple of N c , the dominant gap corresponds to multiple embeddings of the 
Nf = N c order parameter. For Nf = 6, we have 

A^A e *(e ijc + e N(j -3)c). (14) 

This order parameter corresponds to the symmetry breaking pattern 

SU(6) L x SU(6) R x U(l) v -> SU(3) V x U(l) v x U(1) A . (15) 

The £77(3) symmetry is a double embedding of the SU(3) that appear in the Nf = 3 color- 
flavor locked phase. The original U(l)v and approximate U(1)a symmetries are broken, but 
a new U(l)y x U(1)a appears. This symmetry is a subgroup of the original SU (6)l x SU {Q)r 
flavor symmetry which rotates the two 3x3 blocks in equation (|14j). The extra U(l) 
symmetries may be broken by higher order condensates. 

In this section we have analyzed the color-flavor structure of the superconducting ground 
state using a BCS-like energy functional. Another method, which has proven to be very 
useful in the context of liquid 3 He and other systems lly] , is the Landau- Ginzburg functional. 

tin 

This method was applied to color superconductivity in |8J, H2I|. In this case we construct the 
most general energy functional that is consistent with the symmetries of QCD and a simple 
polynomial in the order parameter. In four dimensions, it is usually sufficient to keep terms 
that are at most quartic in the fields. We can analyze the possible ground states by studying 
the minima of the Landau- Ginzburg functional. 

There is one restriction that one has to keep in mind. At T = the free energy of the 
system in not an analytic function of the gap so that, strictly speaking, the free energy 
cannot be expanded as a power series in the order parameter. Only in the vicinity of the 
finite temperature phase transition does the expansion in powers of the order parameter 
have a firm foundation. 

Nevertheless, the Landau-Ginzburg description is very useful in describing a wealth of 
phenomena, even at T = 0. If we restrict ourselves to color and flavor anti-symmetric order 
parameters we can represent the order parameter matrix A?* by the field <f>f, where a is an 
index in the anti-symmetric [N C (N C — l)/2] of color and % in the [Nf(Nf — l)/2] of flavor. 
Color and flavor invariance imply that the Landau-Ginzburg effective potential has the form 

V = -m 2 tr(0 f 0) + Aiftr(0V)l 2 + A 2 tr[(0V) 2 l, (16) 
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where (0t0) ab = (0°)*0*, and m 2 ,Ai and A2 are arbitrary parameters. Of course, we can 
equally well write (|T^|) in terms of (<f><j)')ij = 0"(0")*. Note that there are no det(0) terms, 
because such a term would violate the U(l) of baryon number. The effective potential only 
depends on the eigenvalues of 0*0. This means that the effective potential only depends on 
a number of parameters, N C (N C — l)/2 (or Nf(Nf — l)/2, whichever is smaller). 

For Nf = 2 there is only one quartic term and the minimum occurs for 0° = A5 a3 , 
as we would expect. In the case Nf = 3 there are two possible minima, 0" = A5" and 
0" = A5 a3 5i3. These minima correspond to the color-flavor locked phase and the Nf = 2 
phase. The groundstate depends on the values of the coupling constants. If A2 > 0, the 
ground state is the color-flavor locked phase. Extending the mean field description described 
in the last section to T 7^ 0, we can check that indeed Ai = and A2 > 0, see also flQ. 

For Nf > 3 the effective potential does not depend on the number of flavors. This means 
that there are a large number of ground states all of which are degenerate with an embedding 
of the Nf = 3 color-flavor locked phase. The ground states of the BCS functional discussed 
above are minima of the Landau-Ginzburg potential, but the effective potential (fTHjl does 
not distinguish partially gapped from fully gapped states. In order to construct an effective 
energy functional for which the fully gapped state is the true minimum one has to include 
higher powers in the order parameter. This, of course, would also introduce additional 
parameters and we will not pursue this problem here. 



IV. MORE ON COLOR-FLAVOR-LOCKING IN QCD WITH N c = N f = 3 

In this section we wish to examine the color-flavor-locked state in QCD with three colors 
and flavors in somewhat more detail. For this purpose, we will concentrate on the regime of 
very large chemical potential, \i ^> Aqcd, in which perturbative calculations are possible, 
'he determination of the gap in perturbative QCD has recently attracted some attention 
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22j . The main conclusion is that the gap is dominated by almost collinear 



magnetic gluon exchanges. The magnitude of the gap is 

A ~ 256vr 4 (2/iV / ) 5 /2^- 5 exp ("^) • (17) 

We should emphasize that, strictly speaking, this result contains only an estimate of the 
preexponent. This estimate is obtained by collecting the leading logarithms from electric 
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and magnetic gluon exchanges 18]. There are corrections of order 0(1) that originate from 
improved matching of the gap function at po ~ Ao and po ~ gfi, a self-consistent treatment 
of the Meissner effect, and, possibly, higher order perturbative corrections. Nevertheless, the 
main point is that (|17|) is the result of a well-defined calculation that can be systematically 
improved. 

In previous works, the gap equation was always studied for a one-component color-flavor 
anti-symmetric order parameter. This is appropriate for QCD with two flavors, but in the 
case of more than two flavors the gap equation is more complicated. In the previous section 
we studied a more general N c Nf x N c Nf dimensional gap matrix but restricted ourselves to 
short range, color and flavor antisymmetric interactions. The restriction to antisymmetric 
gap matrices is probably justified for Nf ^ N c , but for Nf = N c the symmetric and anti- 
symmetric order parameters have the same global symmetry, so there is no symmetry reason 
for the symmetric gap parameter to be zero. 

In three flavor QCD, the order parameter has the form 

A# = A A {8«8] - Sffl + A 5 ($?<$* + Sffl. (18) 

We can now repeat the perturbative calculation of the gap using this particular ansatz for 
the order parameter. We will follow the method described in |l8||. Just as in the Nf = 2 
case, the gap depends on frequency and the Dirac structure of the gap matrix is proportional 
to 075(1 + a ■ p)/2. The only new ingredient is that we have to take into account the more 
complicated color-flavor structure when we calculate the Nambu-Gorkov propagator. As in 
the last section, this is most easily accomplished by viewing (fTHJ) as a matrix in a N c Nf 
dimensional color-flavor space. The eigenvalues of this matrix are 

A 8 = A A -A S , A 1 = 2A A + 4A 5 , (19) 

where the subscript indicates the degeneracy. For three degenerate flavors the normal com- 
ponents of the inverse Nambu-Gorkov propagator are proportional to the unit matrix in 
color-flavor space, so they remain diagonal as the anomalous components are diagonalized. 
We can now determine the propagator by inverting the Nambu-Gorkov and Dirac structure 
as in the Nf = 2 case. The off-diagonal propagator S 2 i needed in the gap equation is a 
diagonal matrix in color-flavor space with entries 
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with Aj = (A 8 , A 8 , . . . , Ai) as in (JTHJ). Having determined the propagator we have to 
calculate the color factor. It is given by 



oa = - A {\ A y ab m - ^])(\ a u = -^j^m - s w) (2i) 

for the color-flavor anti-symmetric gap and 

cs = \(\ A f ab (<5?<5j + S^) (X A ) cd = ^1 ($6* + 6?6°) (22) 

for the symmetric gap. The anti-symmetric color factor c A agrees with the result for the 
Nf = 2 order parameter. We can now project the gap equation on the color-flavor anti- 
symmetric and symmetric structures. We find 

9 2 f j i { bfi \ J 2 A A (q ) - A s (q 



(23) 



A MPo) = o / d 9olog i 1 < 5 , = = 

1 2A A (g )+4A s (g ) 



6^ + (2A A (g )+4A 5 (g ))2 
£ 2 f , , / ^ \ J 1 A^(go) - A 5 (g ) 



A 5 (p ) = ^2" y d ?olog 
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JP0"9bl/ I 6 ^q 2 o + (A A (q )-A s (q )) 
1 2A A (g )+4A 5 (g ) 



12 . L? 



g 2 + (2A A (g )+4A s (g ))2 



(24) 



where b = 256tt 4 (2/A^ / ) 5 / 2 ^- 5 . This is a complicated system of coupled equations, but the 
situation simplifies in the weak coupling limit. In this case, we can assume that As <C A A . 
The equation for A A then becomes 

a / \ 9 2 [, , / ¥ \ [2 A A (g ) 1 A A (g ) | 

Aa Po = 2 / log | 1 < - I = = + - , = > • 25 

If it were not for the factor 2 in the denominator of the second term in the curly brackets, 
this would be exactly the same equation as the one we found for the simple Nf = 2 order 
parameter. We can take the factor 2 into account in an approximate way by rescaling the 
integration variable in the second term. In this way we can reduce equation (|23|) to the gap 
equation for the Nf = 2 order parameter with the coefficient b rescaled by a factor 2 -1 / 3 . 
This means that 

A A ~ 2- 1 / 3 256vr 4 (2/iV / ) 5 /V£~ 5 exp ("^J ■ ( 26 ) 
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The equation for As can be analyzed in a similar fashion. We find 



Ac 



<?V21og(2) 



A 4 . 



(27) 



7T 36 

This implies that, formally, As is suppressed by one power of the coupling constant, g. 
In addition to that, we note that the numerical coefficient in (|27)l is quite small, so that 
Ag <C A a even if g is not small. 

These results are easily generalized to an arbitrary number of flavors and colors with 
Nf = N c = N. The eigenvalues of the gap matrix are 



V/V 2 -l 



A A - Ac 



A 1 = (iV-l)A A + (iV+l)A 5 , 



(28) 



and the N c dependence of the anti-symmetric and symmetric color factor is given in (|21|22jh 
From these results we find that the color-flavor-locked gap is given by 



A.. 



2567r 4 (iV - l)" i ^ 1 (2/A^) 5/2 /i^" 5 exp 



7T~ 



6iV 



9 ViV 



Ac = 



9 (N 



( 6N 



1/2 



log(iV - 1) A A . 



(29) 
(30) 



6tt V 2N J \N + I, 

The origin of the various factors of iV in (|2*9*|) is easily explained. The factor in the exponent 
is the color factor that comes from the tree level scattering amplitude of two quarks in a 
color anti-symmetric state. The factor N~ 5 ^ 2 originates from the flavor dependence of the 
screening mass, and the factor (N — 1) raised to the power — (N — 1)/(2N) comes from the 
structure of the color-flavor locked state. This factor implies that for large N, the CFL gap 
is suppressed by a factor with respect to the gap in the Nf = 2 phase. As a consequence, 
for N > 3 the color-flavor locked state (j!8|) is not necessarily the energetically favored state. 
If, on the other hand, we take the large N c limit with Nf = 3 fixed, we find multiple 
embeddings of the iV c = Nf = 3 state. If the large iV c limit is taken with the conventional 
scaling g 2 N c = const, the superconducting gap is suppressed by exp(— y/NZ). This means 
that for very large iV c and Nf fixed, the superconducting ground state is disfavored compared 



to a chiral density wave 
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24|. 



In addition to the gap equation, we would also like to study the condensation energy. To 



order g 2 the grand potential can be calculated from 



8. 25 



n 



{-tr [S(g)£(g)]+trlog [S 1 (q)S(q)}} 



(31) 
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where S(q) and are the Nambu-Gorkov propagator and proper self energy given in 
The grand potential has the form Q = 8/(A§) + /(Ai) where Ag^ are the singlet and octet 
gaps. The functional /(A) is given by 

f(A) = ff 2 [d Po \- A(Po) + ^g + AW-J. (32) 



This result is very similar to the result in the case of short range interactions [8|]. The only 
difference is that the energy dependence of the gap function acts as an explicit cutoff in the 
integrals. We can calculate the integrals using the approximate solution of the Eliashberg 
equation derived by Son Q]. Numerically we find that the condensation energy scales to 
very good accuracy as 

m = (33) 

where A = A(p = 0). This is very similar to the result we found for short range inter- 
actions, equation (jSJ. This suggests that the results obtained in section ITTT1 remain valid in 
the more general case of long range interactions. In particular, we can calculate the energy 
gain of the color-flavor locked state over the Nf = 2 state. The result is only very weakly 
dependent on the gap in a very wide range of A//i. In the weak coupling limit, we find 
e(CFL)/e(Nf = 2) ~ 1.9. This ratio is somewhat smaller than what one would expect based 
on the number of gaps, 9/4 ~ 2.2. 

V. CHIRAL SYMMETRY BREAKING 

The most interesting aspect of the color-flavor locked phase is that chiral symmetry is 
broken, and that the form of the corresponding low energy effective action agrees with QCD 
at low density jj, 2f| • But while the coefficients of the effective lagrangian are complicated, 
non-perturbative quantities in QCD at low density, they can be calculated perturbatively 
at high density. 

In this section we would like to begin this program by calculating the chiral order pa- 
rameter in the color-flavor locked phase. As a first step, we have to calculate the superfluid 
condensate. For a single fermion species we have 

W-Tskl + 5 • m = £ I dpo , 2 A(P °] ■ (34) 
2 71 J -/^ + A(p ) 2 
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In the weak coupling limit, the integrand can be written as a total derivative using the 
differential equation for A(p ) P- 127 1 . The result is 

In the CFL state, the color-flavor structure of the condensate is more complicated. Just like 
the gap matrix, the condensate can be written as 

(C<?75~(1 + a ■ p)^) = MW- ~ + + (36) 

The two condensates 4>a,s can be determined from the octet and singlet gap parameters. 
Not surprisingly, we find that (f>s is small in the weak coupling limit. In the case of (ft a, all 
color and flavor factors drop out and 

The chiral structure of the superfluid order parameter is ipCj^tp = ipRipn — iPl^l- This 
means that pairing takes place between quarks of the same chirality. A convolution of two 
superfluid order parameters (ipRipRji^L^L) will then directly yield the gauge invariant order 
parameter ({iI>l' 1 Pr)(' 1 1>l' 1 1>r)) ■ in the weak coupling limit, factorization is valid and we find 

((Mr)(Mr)) = \ (12^ + 240|) ~ 3<P 2 A . (38) 
In the same way we can calculate many other four fermion operators, like 

<(V^) 2 > = ((Mr)(Mr)) + ($#J>l)$i*J>l)) = H\- (39) 

Here, we have used ((4>liPr)(iPriPl)) — 0. This is a consequence of the chiral structure of 
the superfluid order parameter mentioned above. 

There is an important point that we need to emphasize here. The vacuum expectation 
value ((ijjijj) 2 ) is not an order parameter for chiral symmetry breaking. This is most obvious 
in the case of two flavors. In this case, (rfip, ^ij^rif}) transforms as a 4-vector under the 
chiral SU(2)l x SU(2)r = 0(4). Chiral symmetry restoration then implies 3{(iJjiI}) 2 ) = 
(( , 0ir75'0) 2 ), not ((iptp) 2 ) = 0. In the case of Nj = 3, a true order parameter for chiral 
symmetry breaking is given by 



0i = (^lI^MrI^r), 
14 



(40) 



where A a is a flavor generator. This operator transforms as (8, 8) under SU(3)l x SU(3)r, so 
a non-zero expectation value of Oi definitely implies that chiral symmetry is broken. Fierz 
rearranging 0\, we get an alternative order parameter 



2(N 2 f -1) - - 1 
2 = f N2 L (Mr)(^l) ~ ^WlWrWrWl)- (41) 

f f 



We can calculate the vacuum expectation values of 0\ j2 m the mean field approximation. 
Both turn out to be zero. This does not, of course, imply that chiral symmetry is unbroken. 
The fact that (^1,2) vanishes is due to an accidental symmetry of the mean field approxima- 
tion. The superfluid order parameter is invariant under {Z 2 )l x (Z 2 )r symmetries that act 
on the fermion and anti-fermion fields separately. Because do not have this symmetry, 
they cannot acquire an expectation value in the mean field approximation. This [Z 2 x Z 2 ) 2 
is not a symmetry of QCD. As a result, we expect that 0^ 2 will develop an expectation 
value once higher order perturbative corrections are included. Another way to look at this 
phenomenon is the observation that the color-flavor locked order parameter couples left and 
right handed fields only indirectly, through the vector-like character of the gauge symme- 
try. A non- vanishing expectation value for an operator of the form (LL) (RR) only arises at 
higher order in perturbation theory. 

Since chiral symmetry is broken, we also expect that the standard chiral order parameter 
(■0V) acquires an expectation value. It was noted in ^, 22 1 that the color-flavor locked phase 
has an approximate {Z 2 )l x {Z 2 )r symmetry. If this symmetry were exact, then the quark 
condensate would be zero. In QCD instantons break {Z 2 )l x {Z 2 )r to (Z 2 )v and lead to a 
non-vanishing quark condensate. In 0] we calculated the quark condensate at moderate 
densities, assuming that instantons dominate not only the quark condensate, but also the 
superfluid gap. At very large density, instantons are suppressed and the gap equation is 
dominated by perturbative effects. The results of [2^ are easily generalized to this case. 

In QCD with three flavors, the instanton induced interaction between quarks is given by 
HQS 

C = G 6N C (N 2 - i) e ^3 6 ™3 ( 1^ I \ $L,h *pR,gi ) fefe ^Rm ) ($L,f 3 ipR, g3 ) (42) 

+ 8(A^ 3 + 2) ^ L 'h^ R ^)^LJ 2 <y^R,g 2 )^LJ 3 <y^R,g,) + (L <-> R)\ . 

The coupling constant G is determined by a perturbative calculation of small fluctuations 
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around the classical instanton solution. The result is 



G= dpn(p,n) (2npY 



(43) 



with 



2N C 



Cn 



C 



N 



9' 



p exp 



exp 



8tt 2 

9(p) 2 

0.466 exp(-1.679A^ c )1.34 7V -f 



(44) 



(JV c -l)!(iV c -2)! 

At zero density, the p integral in ()43|) is divergent at large p. This is the famous infrared 
problem of the semi-classical approximation in QCD. At large chemical potential, however, 
everything is under control and G is reliably determined. We find 

5+6 /, _\ 

(45) 



(2k) 



2 6N C (N 2 - 1 



-C N S% Nc A- 5 N 



.5±b ( A 



/ 



ji) \ 2 

Here, A is the QCD scale parameter, b = yA c — |A/ = 9 is the first coefficient of the 
QCD beta function and So = 8n 2 /g 2 . The result shows that, asymptotically, the coupling 
constant G has a very strong power-law suppression ~ /i~( 5+b ) = pT x ^. 

Since G is so small, we can treat the effect of instantons as a perturbation. In the 
color-flavor locked phase the instanton vertex induces a fermion mass term. This can be 
seen by saturating four of the external legs of the interaction (J42|) with the superfluid order 
parameter (J3EJ). Using the results of we find 

6 



(46) 



We note that there are two kinds of fermion mass terms. Both the color singlet structure 
(S ab 5ij) and the color octet structure (5 a i 5 b j) are compatible with the residual SU(3)v sym- 
metry of the color-flavor locked phase. From the result (|46j). we can directly read off the 
singlet and octet mass terms 



m 
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y 



G<P 2 A , 



m 8 



G<Pl 



(47) 



From these results we can also determine the quark condensate. For this purpose we have 
to know the momentum dependence of the fermion propagator in the vicinity of the Fermi 
surface. In principle, this is determined by the momentum dependence of the \x ^ fermion 
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zero mode of the instanton, see 2^| . For simplicity, we assume that the quark mass has the 
same momentum dependence as the gap. In this case, we find 



<« = -2(|l)^H G «. (48) 



From the instanton calculation we can also see how to construct a chiral order parameter 
that is non-vanishing already on the level of the mean field approximation. Consider the 
8-fermion operator 

£> 8 = det(^fl)^*, (49) 
where det^x,^) is a short hand expression for the flavor structure of the instanton vertex 
(142}) . The determinant is invariant under SU(3) x SU(3)r, so Og transforms like the quark 
condensate. From the discussion above it is clear that Og acquires an expectation value in 
the color-flavor locked phase. We find 

(0 8 ) = 60*. (50) 

There is nothing special about Og, other eight-fermion operators are equally good order 
parameters. 

Finally, it is interesting to obtain a few numerical estimates. For definiteness, we will 
consider the chemical potential /i = 500 MeV. Also, we will use A.qcd — 200 MeV. Following 
[19I ] we will be optimistic and use g = 4.2, which corresponds to the maximum of the 
gap as a function of g. In this case we get a substantial gap in QCD with two flavors, 
A (N f = 2) = 130 MeV. The gap in the N f = 2 phase of QCD with three flavors is 
A (iV/ = 3) = 47 MeV. This large reduction comes from the factor NJ 5 ^ 2 in the perturbative 
expression for Ao. This factor is a reflection of the Nf dependence of the screening mass. 
The color anti-symmetric and symmetric order parameters in the color-flavor locked phase 
are A^ = 38 MeV and As = 1.3 MeV. The condensation energy is e = — 32MeV/fm 3 . The 
superfluid order parameter is given by <$a = (144 MeV) 3 . From this we find that the chiral 
condensate is (ipip) = —(35 MeV) 3 , while the expectation value of the four fermion operator 
is significantly bigger {{ipip) 2 ) = (194 MeV) 6 . 



VI. CONCLUSIONS 



In this work we studied the structure of QCD with an arbitrary number of flavors at 
high baryon density. We assumed that the dominant instability of the quark fermi liquid is 
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towards the formation of a pair condensate. In order to study the ground state, we introduced 
a renormalized mean field grand canonical potential which determines the condensation 
energy as a function of the color-flavor structure of the gap matrix. This potential is not 
derived from QCD, but we expect it to correctly represent the essential dynamics of the 
high density phase. From the analysis of the grand potential we find that for any number 
of flavors greater than two, chiral symmetry is broken while a vector-like flavor symmetry 
remains. The most interesting case is QCD with Nf = 3 flavors. In this case, the symmetry 
breaking pattern at high density matches the one at low density. The color-flavor locked 
state in Nf — 3 is also distinguished by the fact that it has the biggest condensation energy 
per flavor, and that it leaves the largest subset of the original flavor symmetry intact. 

In the second part of this work we studied the Nf = 3 phase in more detail. We calcu- 
lated the color symmetric and anti-symmetric order parameters in perturbative QCD. We 
showed that A s is suppressed by one power of g in the weak coupling limit, and that the 
color-flavor locked state is the energetically favored state. We calculated (ipip) and ((4>tp) 2 ). 
We explicitly showed that the chiral condensate is suppressed because of an approximate 
Z2 x Z2 symmetry. In general, the expectation value of four-fermion operators is not sup- 
pressed, but all gauge invariant four-fermion chiral order parameters vanish in the mean 
field approximation. Non- vanishing gauge invariant chiral order parameters can be obtained 
by considering higher dimension operators, semi-classical (instanton) effects, or higher order 
perturbative corrections. 
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VII. ERRATUM 

There is a mistake in equs. (11) and (12) in which we describe the symmetry breaking 
pattern in the case N c = 3, Nf = 4. The second set of 577(2) generators given in equ. (12) 
does not generate a symmetry of the order parameter (and it does not commute with the 
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first set of generators). The correct generators are 



(51) 



Note that the second SU (2) is a pure flavor symmetry. As a consequence, the symmetry 
breaking pattern is 



SU{A) L x SU(A) R x U{l) v -> SU(2) V x SU(2) V x SU(2) A 



(52) 



and the axial SU {A) a symmetry is not completely broken. I thank H. Murayama for pointing 
out this error. 

There is a factor 2 mistake in equation (31) for the thermodynamic potential. To o(g 2 ) 
the thermodynamic potential is given by 



Q = ^/7Sl{- tr [^) S ^)]+trlog[5 - 1 (g)5(g)]} 



(2nY 
1 r d 4 q d 4 k 
A J (2tt) 4 (2tt) 41 



s( q )r«s( q + k)r b \ D;l(k), 



(53) 
(54) 



where is the quark gluon vertex function. The last term can be simplified using the gap 
equation. We find 

1 f d A q f 1 



{} = 27 (2^i~2 tr[%)S(g)]+trl ° g ^ 1(g)5(g) 



(55) 



Note that the first term differs from equ. (31) by a factor 1/2. As a consequence, equ. (32) 
should read 



/(A) = 



A(po) 2 



2^ + A(po): 



+ a/^ + A^o) 2 - Pol- 



and equ. (33) is replaced by 



/(A) 



47T 2 ° 



(56) 



(57) 



which is identical to the result in BCS theory. 

Finally, equ. (48) gives an estimate of the quark condensate in the CFL phase. The 
complete weak coupling result is 

,2 \ 



(W) = "2 



2n 2 



12G01, 



(58) 



which was obtained in T. Schafer, Phys. Rev. D65, 094033 (2002). 
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N c 


N, 


Npar 


gaps (deg) 


A 


-e/(N c N f ) 




rank 


3 


2 


3 


A (4), (2) 


A 


eo 


6 


2 


3 


3 


9 


A (8), 2A (1) 


0.80A 


1.27eo 


8 


2 


3 


4 


18 


A (8), 2A (4) 


0.63A 


1.21eo 


6 


2 


3 


5 


30 


A (5), 2A (7), 3A (3) 


0.43A 


1.18eo 


3 


1 


3 


6 


45 


A (16), 2A (2) 


O.8OA0 


1.27eo 


9 


3 



TABLE I: Spectrum and symmetry properties of the s-wave superfluid state in QCD with N c = 3 
colors and Nf flavors. N par = N C (N C — l)Nf(Nf — 1)/4 is the number of totally anti-symmetric gap 
parameters. The fourth column gives the relative magnitude of the gaps in the fermion spectrum, 
together with their degeneracy. The next two columns give the magnitude of the gap and the 
condensation energy per species in units of the Nf = 2 values. These ratios are independent of the 
coupling in the weak-coupling limit. The last two columns show the dimension and the the rank 
of the residual symmetry group. 
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